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ABSTRACT 

Aims. In the effort to understand the link between the structure of galaxy clusters and their galaxy populations, we focus on MACSJ1206.2-0847 
at z ~ 0.44 and probe its substructure in the projected phase space through the spectrophotometric properties of a large number of galaxies from 
the CLASH-VLT survey. 

Methods. Our analysis is mainly based on an extensive spectroscopic dataset of 445 member galaxies, mostly acquired with VIMOS@VLT as 
part of our ESO Large Programme, sampling the cluster out to a radius ~ 2/?200 (4 Mpc). We classify 412 galaxies as passive, with strong Hd 
absorption (red and blue galaxies), and with emission lines from weak to very strong. A number of tests for substructure detection are applied to 
analyze the galaxy distribution in the velocity space, in 2D space, and in 3D projected phase-space. 

Results. Studied in its entirety, the cluster appears as a large-scale relaxed system with a few secondary, minor overdensities in 2D distribution. 
We detect no velocity gradients or evidence of deviations in local mean velocities. The main feature is the WNW-ESE elongation. The analysis of 
galaxy populations per spectral class highlights a more complex scenario. The passive galaxies and red strong Hd galaxies trace the cluster center 
and the WNW-ESE elongated structure. The red strong Hd galaxies also mark a secondary, dense peak ~ 2 hj^ Mpc at ESE. The emission line 
galaxies cluster in several loose structures, mostly outside ^ 200 - Two of these structures are also detected through our 3D analysis. The observational 
scenario agrees with MACS J1206.2-0847 having WNW-ESE as the direction of the main cluster accretion, traced by passive galaxies and red 
strong Hd galaxies. The red strong Hd galaxies, interpreted as poststarburst galaxies, date a likely important event 1-2 Gyr before the epoch of 
observation. The emission line galaxies trace a secondary, ongoing infall where groups are accreted along several directions. 

Key words, galaxies: clusters: individual: MACS J1206.2-0847 - galaxies: clusters: general - galaxies: kinematics and dynamics - galaxies: 
evolution - cosmology: observations 


1. Introduction 

In the hierarchical scenario for large-scale structure formation, 
clusters of galaxies are not relaxed structures. Numerical simu¬ 
lations show that clusters form preferentially through anisotropic 
accretion of subclusters along the large-scale structure (LSS) fil¬ 
aments (Colberg et al. ll998lll999l and references therein). From 
the observational point of view, it is well known that a large frac¬ 
tion of clusters (30%-70%) contain substructures, as shown by 
studies based on optical (e.g., Baier & Ziener W911\ Geller & 
Beers [T9821 Girai'di et al. [T9^ Ramella et al. l2007l Wen & 
Han l2013ll . on X-ray data (e.g., Jones & Forman [19991 Buote 
120021 Jeltema et al. 120051 Zhang et al. 120091 Bohringer et al. 
12010b . and on gravitational lensing techniques (e.g., Athreya 
et al. I 2 OO 2 I Dahle et al. l2002l Smith et al. [2005] Grillo et al. 
I2014ab . indicating that past signatures of cluster accretion are 
quite common. Much progress has been made in the study of 
cluster accretion phenomena (see Feretti et al. 120021 for a gen¬ 
eral review) and recent, dedicated studies have often focused on 
a few major, ongoing cluster mergers, for instance in the con- 
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text of Dynamical Analysis Radio Clusters project (DARC, e.g., 
Girardi et al. l201 lb and MUlti-Wavelength Sample of Interacting 
Clusters project (MUSIC, Maurogordato et al. l201 lb . Other ded¬ 
icated studies have focused on larger scales, i.e., on cluster ac¬ 
cretion through filaments (e.g., Fadda et al. 120081 Dietrich et al. 
12012b . On-going cluster formation is also evident in distant clus¬ 
ters as pioneering studies have shown that most clusters identi¬ 
fied at z > 0.8 are elongated, clumpy, or with a filamentary struc¬ 
ture (e.g., Donahue et al. 119981 Gioia et al. 1 19991 Demarco et al. 
120071 Fassbender et al. l201 lb . 

It is well established that the properties of cluster galaxies 
differ from those of field galaxies and that clusters are character¬ 
ized by radial gradients. Galaxies in denser, central regions are 
usually of earlier morphological type, redder color, and lower 
star formation rate (hereafter SFR; e.g., Gerken et al. 12004b . 
However, the precise details of the connection between galaxy 
evolution and cluster environment are still unknown. Several 
physical mechanisms can be involved in modifying galaxies in 
the cluster environment (see Fig. 10 by Treu et al. l2003l and refs, 
therein). Moreover, it is clear that a scenario in which galaxies 
are accreted in a static cluster environment is too simplistic. In 
fact, there is evidence of a connection between the fraction of 
active galaxies and the presence and/or the amount of substruc- 
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ture, generally interpreted as an enhanced star formation (SF) in 
member galaxies caused by the process of cluster mergers (e.g., 
Caldwell & Rose 1997; Ferrari et al. l2005l Owen et al. f2005]) . 
An alternative interpretation is that quenching SF is less effec¬ 
tive in just forming multicomponent clusters (see Cohen et al. 
120141 for further discussions). From a theoretical point of view, 
during a major cluster merger, the SF in gas-rich galaxies can 
decrease or, alternatively, be triggered depending on the specific 
interaction between the intracluster medium (ICM) and the in¬ 
terstellar gas within the galaxies (e.g., Fujita et al. 119991 Bekki 
et al. 120101) . Time-dependent tidal fields of merging groups and 
clusters of galaxies can also trigger secondary starbursts in their 
member galaxies (Bekki IT999t . 

In this context, poststarburst galaxies (PSBs, Dressier & 
Gunn 1198^ play a crucial role. In fact, their spectral signature, 
i.e., strong Balmer absorption lines and weak or no emission 
lines, is generally attributed to a vigorous SF that has recently 
decreased significantly or ceased altogether. More properly, one 
can think of a real poststarburst phase or, alternatively, a star¬ 
forming truncated phase (e.g.. Dressier et al. 120091) . The last SF 
event in PSBs galaxies ended a few Myr to 2 Gyr before the 
time of observation, on the basis of the lifetime of the stars re¬ 
sponsible for the strong Balmer lines (e.g., Poggianti & Barbaro 
[19971 Balogh et al. [19991 Poggianti et al. [T999l Mercurio 
et al. 120041 hereafter M04). More specifically, the typical time 
elapsed since the last SF in the red PSBs lies in the range 1- 
2 Gyr, while in the blue PSBs the interruption of SF occurred 
within the last 0.5 Gyr (M04). The location of PSBs in the par¬ 
ent cluster can give important information on the driving phys¬ 
ical mechanism. For instance, the synchronized triggered SF in 
the cluster merging simulations of Bekki et al. (120101 ) gives ori¬ 
gin to a population of PSBs having a spatial distribution differ¬ 
ent from the rest of the galaxy population. Analyzing nine clus¬ 
ters at z ~ 1 drawn from the GCLASS survey, Muzzin et al. 
(12014b have found that PSBs appear to trace a coherent ring in 
the projected phase space. Alternatively, Dressier et al. (12013b 
have found that in the IMACS Cluster Building Survey (five rich 
clusters at 0.31 < z < 0.55) the PSB population traces the same 
spatial distribution of the passive population. 

Recently, significant interest has focused on cluster outskirts, 
too. In fact, although there is evidence that SF is suppressed 
on average outside 2 R 20 cEl in both low- and high-redshift clus¬ 
ters (e.g., Lewis et al. 120021 Gomez et al. 120031 Rasmussen et 
al. 12012b . several studies have associated starburst features with 
infalling groups (e.g., Owen et al. 120051 M04; Poggianti et al. 
120041 Oemler et al. 12009b . In nearby clusters an increase of SF 
along the filaments feeding the clusters has been found (e.g., 
Fadda et al. 12008b . This seems particularly true around 1.5-2^200 
and for galaxies belonging to groups, suggesting that a rela¬ 
tively high density in the infalling regions promotes interactions 
between galaxies (e.g., galaxy-galaxy harassment) and follow¬ 
ing momentary bursts of SF (Scott et al. 120071 Mahajan et al. 
12012b . Recent simulations suggest that cluster galaxies may be 
pre-processed before the infall into the cluster environment, al¬ 
though the postprocessing is important as well (Vijayaraghavan 
& Ricker l2nT3b . 

Optical multiband data and particularly multiobject spec¬ 
troscopy are a consolidated way to investigate cluster substruc¬ 
ture and cluster merging phenomena (e.g., Girardi & Biviano 
12002b . The few recent studies of clusters with hundreds of spec¬ 
troscopic members have shown the power of using large spectro- 


* The radius Rg is the radius of a sphere with mass overdensity S 
times the critical density at the redshift of the galaxy system. 


scopic catalogs to study cluster internal structure (e.g., Owers et 
al. l20111 Munari et al. 12014b . In some studies, the correlation of 
passive and active galaxies populations with the substructure of 
the parent cluster has been investigated (e.g., Czoske et al. 120021 
M04; Mercurio et al. 120081 Oemler et al. 120091 Ma et al. 12010b . 
Very large spectroscopic datasets are needed to study cluster 
substructure in phase space for different spectral types and to 
infer cluster assembly history. The present study represents a pi¬ 
lot study in this new direction. It focuses on the galaxy clus¬ 
ter MACS J1206.2-0847 (hereafter MACS 1206; Ebeling et al. 
12009b at z ~ 0.44, discovered and first described in the REFLEX 
Galaxy Cluster Survey catalog (RXC J1206.2-0848; Bohringer 
et al. 12004b . 

1.1. MACS1206 and the CLASH-VLTprogramme 

The subject of this study, the galaxy cluster MACS 1206 is part 
of the survey of 25 target clusters of the HST Multi-Cycle 
Treasury program “Cluster Lensing And Supernova” (CLASH, 
RL; M. Postman; Postman et al. 12012b . Our study is based 
on an extensive spectroscopic dataset collected within the ESO 
Large Programme 186.A-0798 “Dark Matter Mass Distributions 
of Hubble Treasury Clusters and the Eoundations of ACDM 
Structure Eormation Models” (hereafter CLASH-VLT, PL: P. 
Rosati; Rosati et al. 12014b . This program aims at obtaining spec¬ 
tra for at least 500 cluster members for each of the 14 southern 
CLASH clusters with 0.2 < z < 0.6, using the Visible Multi- 
Object Spectrograph (VIMOS) at VLT (CLASH-VLT survey). 
The study of the combined evolution of galaxies and parent clus¬ 
ter requires i) the spectroscopic identification of a large number 
of cluster galaxies in each cluster; ii) sampling a region out to 
at least twice R 200 ', and iii) a range of cluster redshifts covering 
a relevant epoch cluster formation. Specifically, the past 5 Gyr 
play a crucial role in cluster assembly history since the over¬ 
all morphological content, the fraction of galaxies of different 
spectral types, and the SE history of cluster galaxies are known 
to significantly evolve since z ~ 0.5 (e.g.. Butcher & Oemler 
[T9841 Co uch & Sharpies [T9g7l Dressier et al. [19971 Eritz et 
al. 120051 Poggianti et al. 12009b . Useful complementary infor¬ 
mation from optical. X-ray and radio data, and high-resolution 
HST imaging is available, too. Other CLASH-based studies, fo¬ 
cused on other aspects of MACS 1206, have already been pre¬ 
sented and are referred to and discussed throughout this paper 
(Umetsu et al. 120121 Zitrin et al. 120121 Biviano et al. 120131 
B13; Annunziatella et al. 120141 Grillo et al. 12014bl Presotto et 

The paper is organized as follows. We present our catalog 
and galaxy spectral type classification in Sects. 2 and 3, respec¬ 
tively. Section 4 concentrates on the analysis of cluster structure. 
Section 5 is devoted to the interpretation and discussion of our 
results. We give our summary and conclusions in Sect. 6. 

Unless otherwise stated, we give errors at the 68% con¬ 
fidence level (hereafter c.l.). Throughout this paper, we use 
Hq - 70 km s“' Mpc^' in a flat cosmology with Go = 0.3 and 
Ga = 0.7. In the adopted cosmology, 1' corresponds to 0.341 
hjQ Mpc at the cluster redshift. 

2. Redshift catalog and member selection 

The full redshift catalog of MACS 1206 consists of 3292 ob¬ 
jects, with measured redshifts mostly acquired as part of our 
ESO Large Programme 186.A-0798 (P.I.: Piero Rosati) using 
the VIMOS@VLT/UT3. Additional archival VIMOS data have 
been homogeneously reduced from programs 169.A-0595 (PL: 
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Hans Bohringer) and 082.A-0922 (RI. Mike Lerchster) and have 
provided 952 spectra in the cluster field. All data are reduced 
with the same VIPGI package (see Scodeggio et al. l2005t . Four 
redshift quality classes were defined “secure” (QF=3), “based 
on a single-emission-line” (QF=9), “likely” (QF=2), and “in¬ 
secure” (QF=1), having > 99.99%, ~ 92%, ~ 75%, and < 40% 
probability to be correct, respectively. Further details on the ESO 
Large Programme 186.A-0798 are reported elsewhere [Rosati et 
al. (I2014I I and in Rosati et al. (in prep.)]. The catalog also in¬ 
cludes 22 spectra acquired with FORS2@VLT as part of ESO 
Programme ID 089.A-0879 (PI.; Raphael Gobat) and additional 
literature redshifts taken from Lamareille et al. (120061 1 object), 
Jones et al. (120041 3 objects), Ebeling et al. (120091 25 objects), 
and 21 objects observed with IMACS-GISMO at the Magellan 
telescope (Daniel Kelson private communication). The full cata¬ 
log, with the exception of QP=1 redshifts has already been used 
by B13 and Lemze et al. (120131 1 and is publicly available. In this 
paper, we do not consider redshifts with QP=1 and 2. Our refer¬ 
ence sample includes 1920 galaxies. 

Biviano et al. (I2013I I quantified the completeness of the spec¬ 
troscopic catalog by computing the ratio between the Rc-band 
number counts of objects with measured redshift and all photo¬ 
metric objects. In particular, B13 show that in the virial region 
the completeness is ~ 0.6 at Rc ~ 20, ~ 0.3 at Rq ~ 23, and 
then declines very sharply (see their Pig. 4). However, although 
the redshift catalog covers a fraction of the photometric objects, 
it is essentially unbiased with respect to galaxy color. Biviano 
et al. (I2013I I also evaluated the spatial completeness of the spec¬ 
troscopic sample finding a good uniformity with a mild radially- 
dependent incompleteness. The completeness varies with radius 
from 0.6 at the cluster center to 0.2 at the external cluster limit 
(see their Pig. 5). Only two low-overdensity features are shown 
in their Pig. 5; the first in the center, elongated towards the south 
direction, the second 10' to the east. Our sample and that of 
B13 do not differ in their spatial distributions according to the 
2D Kolmogorov-Smirnov test (hereafter 2DKS-test, Pasano & 
Pranceschini 119871 1. A posteriori, we verified that none of the 
galaxy clumps or subclusters we discuss in this paper resem¬ 
ble the specific overdensity features detected in the complete¬ 
ness map by B13 (cf. their Pig. 5 with our results and figures in 
Sect. 3.2). We conclude that our results are not affected by the 
redshift spatial incompleteness. 

To select cluster members, we applied the two-step method 
called “peak-Hgap” (Ph-G) already applied by B13 and previous 
studies (e.g., Girardi et al. 1201 H and refs, therein). The method 
is a combination of the ID adaptive-kernel method DEDICA 
(Pisani I1993I I and the “shifting gapper”, which uses both posi¬ 
tion and velocity information (see Padda et al. 119961 Girardi et 
al. [T996] |. Por the center of MACS 1206, we adopted the posi¬ 
tion of the brightest cluster galaxy [BCG, R.A.=12^06™12? 15, 
Dec.=-08°48'03.3" (J2000)]. In the first step, we detected 

MACS 1206 as a peak at z ~ 0.4385 populated by 466 galax¬ 
ies. The second step leads to 445 cluster members. 

By applying the biweight estimator to the 445 cluster mem¬ 
bers (Beers et al. 119901 ROSTAT software), we computed a mean 
cluster line-of-sight (LOS) velocity {V} = (cz) = (131 843+49) 
km s“', corresponding to a mean cluster redshift (z) = 0.4398 + 
0.0002. We estimated the LOS velocity dispersion, cry, by us¬ 
ing the biweight estimator and applying the cosmological cor¬ 
rection and the standard correction for velocity errors (Danese 
et al. 119801) . We obtained cry = 1035];45 km s“\ where errors 
are estimated through a bootstrap technique. Both (V) and cry 
are in agreement with the values obtained by B13 using a larger 
galaxy sample. Out of the reference sample of 445 cluster mem¬ 


bers (hereafter the TOT sample), we also considered the sample 
of 249 galaxies within R 200 (hereafter the R200 sample). The es¬ 
timate of R 200 spans the range 1.96-2.08 Mpc, with a typical 
Icr uncertainty of ~ 5%, depending on the method of analysis 
(see Table 2 of B13). Hereafter the value R 200 = 2 Mpc is 
assumed. 

Throughout the paper we also use color and morphological 
information coming from the Subaru Suprime-Cam data. These 
data were retrieved from the SMOKA archival and reduced 
following the same steps as adopted by Nonino et al. (120091 1. 
Purther information can be found in Presotto et al. (120141) and in 
Mercurio et al. (in prep.), where this photometric information is 
more extensively used. 

3. Spectral type classification and galaxy 
populations 

3.1. Spectral type classification 

Out of a total of 445 member galaxies, we have the spectra 
of 415 galaxies acquired at VIMOS@VLT (405 objects) and 
PORS2@VLT (10 objects), other redshifts coming from the lit¬ 
erature or other sources (see Sect.|2]). Poliowing most of the liter¬ 
ature, we used emission lines and H(5 absorption lines to classify 
galaxies. We used wavelength-calibrated and flux-calibrated 
galaxy spectra as obtained from VIPGI (Scodeggio et al. l2005t . 
We corrected for galactic extinction following Schlegel et al. 
(119981) and using the IRAP0task DEREDDEN included in the 
package ONEDSPEC. Pinally, the calibrated and dereddened 
spectra were corrected for the measured velocity dispersion by 
using DISPCOR. 

We measured the equivalent widths (EWs) for the emission 
lines [Oil], [OIII], and when available Her. All the spectra cover 
the [Oil] region with the exception of 18 galaxies. We also mea¬ 
sured the EW for Hd (see M04 for the definition of wavelength 
ranges). We detected 185 galaxies with evidence of emission 
lines. The 185 galaxies with emission lines were divided in four 
classes depending on the strength of the [Oil] emission line, or of 
the [OIII] emission line for the four cases where the [Oil] region 
is not covered. In these four cases we assumed a [OIII]/[OII] 
flux ratio equal to one. We considered the wELG, mPLG, sELG, 
and vsELG classes, i.e., those formed of galaxies with weak, 
medium, strong, and very strong emission lines (see Table [TJ. 
Pigure [1] shows the co-added spectra for emission line galaxies 
and other classes. All the emission line galaxies, with the excep¬ 
tion of wELGs, show similar distributions of the clustercentric 
distances (see Sect. 13.2b . and we aggregated them in the LEG 
class. 

Among the non-emission line galaxies, the 44 galaxies with 
EW(Hd) >3A were classified as strong Hd absorption galaxies 
(HDS sample), interpreted as PSBs or galaxies with truncated 
SP (see Sect.[TJ. Compared to red galaxies, blue galaxies require 
larger EW(Hd) to be identified as galaxies that have a peculiar 
star-forming activity. Thus, a more precise way to classify strong 
Hd galaxies uses the diagram of EW(Hd) vs (B — R) color or, 
alternatively, vs the strength of the 4000 A break, which corre¬ 
lates with color (Couch & Sharpies 119871 Barger et al. 119961 
Balogh et al. 11999b . We considered 38 strong Hd galaxies with 

^ http://smoka.nao.ac.jp/SUPsearch 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Table 1. Spectral type classification. 


Galaxy characteristics 

Class 

Ng 

EW([OII])‘‘ 

EW(H(5) 

Passive 

PAS 

183 

absent 

absent 

Strong Hh absorption & red color 

HDSr*’ 

38 

absent 

> 3A 

Strong H5 absorption & blue color 

HDSb*’ 

5 

absent 

> 5A 

Weak emission lines 

wELG 

17 

> -7A 

any 

Medium emission lines 

mELG" 

24 

(-15A, -7A] 

any 

Strong emission lines 

sELG" 

101 

(-40A,-15A] 

any 

Very strong emission lines 

vsELG" 

43 

< -40A 

any 


Notes, (a) In four cases EW([OIII]) were used, see the text and Table|2j (b) The HDS class (Ng = 44) is formed of HDSr and HDSb galaxies, and 
an additional galaxy with strong absorption, which can be classified as neither HDSr nor HDSb; (c) The ELG class (Ng = 168) is formed of 
mELG, sELG, and vsELG. Throughout the paper we use PASs, HDSrs, HDSbs, and w/m/s/vsELGs to refer to the respective member galaxies. 



4000 4500 5000 5500 

Wavelength [A] 


Fig. 1. Co-added (rest-frame) spectra of cluster galaxies per 
spectral class. The positions of [Oil] and H(5 lines used in the 
classification are indicated. 

red colors and EW(H5)> 3A (HDSr) and 5 strong galaxies 
with blue colors and EW(Hb)> 5A (HDSb). In this separation, 
we followed M04 using their threshold in the differential color 
(B - Rc)dis - -0.5, where the differential color is defined as 
the observed color minus the color of the color-magnitude rela¬ 
tion at the given magnitude. The color-magnitude relation was 
obtained as fitted to the spectroscopically confirmed members 
(B - Rq - 3.587 - 0.0714/?c), where Rq indicates the Kron 
magnitudes and the B - Rc color is computed on aperture mag¬ 
nitudes (Mercurio et al. in prep.). In their original study, M04 
based their HDS classification scheme directly on the strength 
of the 4000 A break, while we use the corresponding color. The 


color-magnitude correction allows us to take into account the 
large magnitude range spanned by our dataset. 

With the exception of three galaxies whose signal-to-noise 
is too low according to our standards (S/N~2 for the Hb absorp¬ 
tion line), the remaining 183 galaxies were classified as passive 
(PAS). We note that the BCG is classified as a passive galaxy 
in spite of the presence of a [Oil] emission line in the VIMOS 
spectrum since the HST photometric data suggest that this line 
can be associated with a blue compact source and/or peculiar 
features blended with the BCG (see Presotto et al. 120141 for the 
relevant discussion). 

Summarizing, we have spectral type classifications for 412 
galaxies, 217 galaxies within f? 200 - The spectroscopic catalog is 
electronically published in Table |2l available at CDS. Table [T] 
summarizes the numbers of galaxies assigned to each spectral 
class. 

To compare our spectral-classification statistics to that of 
Dressier et al. (120131 1. which considers five rich clusters at 
0.31 < z < 0.54, we also computed the fractions per spectral 
class within 1.5 R 200 and Rq < 22.3. We obtained the fraction 
(PAS-hHDS)/TOT= 67% and HDS/PAS=24%. These values do 
not change significantly when taking into account the radial spa¬ 
tial incompleteness (see Pig. 5 of B13 and our Sect.|2]l, i.e., nor¬ 
malizing the number of galaxies in different cluster radial bins 
(R < 0.4 Mpc, 0.4 < /? < 1 Mpc, I < R < 2 Mpc, 
R > 2 Mpc) using different spatial completeness correc¬ 
tions (i.e., dividing by 0.6, 0.5, 0.4, 0.3). Our values are in good 
agreement with the corresponding fractions (PAS-i-PSB)/TOT 
and PSB/PAS reported by Dressier et al. (120131 see the cluster 
values in their Pig. 16 and Table 4). 

We also used some results from the morphological analy¬ 
sis of our galaxy sample as derived from Suprime-Cam@ Subaru 
data (Mercurio et al., in prep.). The PAS, HDSr, wELG, and 
ELG sequences of spectral type populations correspond to bluer 
and bluer galaxies (see Pig.|2]). Por each spectral class with the 
exception of HDSb, Pigure [3] shows the distribution of Sersic 
index n, as determined in the /^cband Subaru image using the 
GALAPAGOS software (Barden et al. l2012l l and a single Sersic 
profile. The value of n is larger for more prominent bulges and 
n-2-2.5 can be considered the transition value between disk- 
dominated and bulge-dominated galaxies (e.g., M04; Barden et 
al. 120051 Pisher & Drorv l2008] l. We note a good correlation be¬ 
tween the Sersic index n and spectral classes, with PASs mostly 
bulge-dominated and ELGs mostly disk-dominated, while HDSr 
and wELG galaxies have intermediate values of n. 
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Table 2. Spectral classification of galaxies in the cluster MACS J1206.2-0847 


ID 

a(J2000) 

(5(12000) 

Z 

B-Rc 

EW([OII]) 

A 

EW(H(5) 

A 

Class 

13889 

19586 

12 : 06 : 55.15 
12 : 06 : 05.16 

-08 : 56 : 29.7 
-08 : 53 : 20.6 

0.4367 

0.4360 

1.77 ±0.01 
1.26 ±0.01 

-1.50 ±0.30 
-2.80 ±0.50 

6.90 ± 0.20 
2.50 ± 0.50 

wELG 

wELG 


Notes. This table will be made available at the CDS. A portion is shown here for guidance regarding its form and content. Column 1: running ID 
for galaxies in the presented sample; Cols. 2 and 3: R.A. and Dec. (J2000); Col. 4: spectroscopic redshift, z\ Col. 5 B - Rc color; Cols. 6-7: [Oil] 
and equivalent widths. We list ‘....’ when the EW cannot be measured because the line does not lie within the available wavelength range. We 
list ‘0.0’ for EW[OII] when no emission is observed. Column 8: the spectral classification. 
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Fig. 2. B - Rc vs Rq color-magnitude diagram with our seven 
spectral classes. The black solid line is the color-magnitude re¬ 
lation, {B - /?c)diff = 0. The blue dashed line indicates the 
{B - f?c)diff = -0.5 value which we use to separate HDSr and 
HDSb following M04. The two red dashed lines indicate the lo¬ 
cus of the red sequence galaxies, |(Z? - /?c)difrl < 0.3. 


3.2. Space and velocity distributions of galaxy populations 

Figure |4] shows that PAS and HDSr galaxies are more spatially 
clustered than ELGs. To make a quantitative comparison among 
different galaxy populations, we applied the Kruskall-Wallis 
test to the clustercentric distances R (KW-test, e.g., Ledermann 
11982b . This test is a non-parametric method for testing whether 
samples originate from the same distribution. The KW test leads 
to significant results (> 99.99% c.l.); at least one of the sam¬ 
ples is different from the other samples. In Fig.|5]there is a clear 
dichotomy between the population formed of PAS, HDSr, and 
wELG galaxies and the population formed of mELG, sELG, 
and vsELG galaxies, the galaxies of the former group being 
more clustered. We also compared spectral classes two by two. 
We applied the ID Kolmogorov-Smirnov test (hereafter IDKS- 
test; e.g., Ledermann I 1 982b to compare the distributions of clus¬ 
tercentric distances. When comparing one of the PAS, HDSr, 
wELG classes with one of the mELG, sELG, vsELG classes 


Fig. 3. Cumulative distribution of Sersic index of galaxies per 
spectral class. 


we detected significant differences (at > 98% c.l.). The com¬ 
parison between PAS and HDSr, PAS and wELG, and HDSr 
and wELG shows no differences, and similarly the comparison 
between mELG and sELG, mELG and vsELG, and sELG and 
vsELG. The results of the IDKS-test confirm the dichotomy be¬ 
tween the population formed of PAS, HDSr, and wELG galax¬ 
ies and that formed of mELG, sELG and vsELG galaxies and 
prompted us to treat mELGs, sELGs, and vsELGs together in 
the combined ELG class, separately from wELGs. This is a more 
detailed view, across several spectral types, of the well-known 
spatial segregation of spiral/elliptical galaxies and of blue/red 
galaxies in local and distant clusters (e.g., Melnick & Sargent 
119771 Dressler ri980l Whitmore et al. 119931 Abraham et al. 119961 
Dressier et al. 1 19991 Gerken et al. 12004b . 

In a second-order analysis, we applied the 2DKS-test to com¬ 
pare the distributions of projected positions. This is particularly 
meaningful for the MACS 1206 dataset, which is large enough to 
search for differences in the internal structure. This test is more 
sensitive than the above IDKS-test and, in addition to confirm¬ 
ing the above dichotomy, it is able to detect differences within 
the ELG class, specifically between mELG and vsELG, and be¬ 
tween sELG and vsELG (marginally, at the < 95% c.l.). The 
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Fig. 4. Spatial distribution of the the 445 cluster members high¬ 
lighting the spatial segregation between early- and late-type 
galaxies across our spectral classes. Each classified galaxy is 
marked by a symbol: PAS (small red circles); HDSr (magenta 
crosses); HDSb (large blue crosses); wELG (green triangles); 
mELG (blue triangles); sELG (blue squares); vsELG (blue cir¬ 
cles). The circle centered on the BCG encloses the R 200 region. 

Table 3. Kinematical properties of the whole cluster and galaxy 
populations. 


Sample 

Ng 

< V > 

km s ' 

(Tv 

km : 

TOT 

445 

131843 ±49 

1035 

R200 

249 

131776 ±74 

1169 

PAS 

183 

131980 ±69 

940 

HDSr 

38 

131630 ±196 

1188 

HDSb 

5 

132836 ± 767 

1307 

wELG 

17 

131502 ±307 

1215 

ELG“ 

168 

131690 ±78 

1014 

mELG 

24 

132121 ±192 

917 

sELG 

101 

131698 ±99 

996 

vsELG 

43 

131485 ±170 

1102 


+27 

45 

+54 


:!I9 

-141 

+415 

i3 


-161 


Notes. Column 1: sample ID [(a) The ELG class is formed of mELG, 
sELG, and vsELG]; Col. 2: the number of assigned galaxies, Ng', Cols. 3 
and 4: mean LOS velocity and velocity dispersion of galaxies, (V) and 
(Tv, with respective errors. 


combined mELG-PsELG class differs from vsELG at the 96% 
C.I.. The observed difference is related to the cluster substructure 
and will be further analyzed and discussed in following sections. 

The application of the KW-test to the rest-frame velocities 
Vrf [i.e., {V - (E))/(l + z)] gives no significant result, i.e., the 
velocity distributions of different spectral classes can originate 
from the same parent distribution (see Eig.|6] upper panel). We 
also applied the KW-test to the absolute value of Erf in such a 
way to be more specifically sensitive to differences in cry val- 


Fig.5. Cumulative distributions of clustercentric distance R of 
galaxies per spectral class, showing the spatial segregation from 
early- through late-type galaxies. 


ues: no significant result was obtained (see Eig.j^ lower panel). 
When comparing two by two the (V) estimates of different spec¬ 
tral classes (see Table O, we found no significant difference. 
The same is true for cry. Our results corroborate previous stud¬ 
ies claiming for a velocity distribution which does not depend 
dramatically on galaxy color or spectral type (Rines et al. 120051 
Mahajan et al. 120111 Rines et al. I2013II . including the results 
of B13 on MACS 1206 itself. On the other hand, several pre¬ 
vious studies have found significant differences in the velocity 
distributions of different galaxy populations, the velocity disper¬ 
sion of the population of blue SE galaxies being higher than that 
of the population of red passive galaxies (e.g., Tammann ll972l 
Moss & Dickens 119771 Sodre et al. 119891 Biviano et al. 119921 
Zabludoff & Eranx|T993] Colless & Dunn lT996l Biviano et al. 
119971 Adami et al. 119981 Dressier et al. I1999I I. The cause for 
this controversy may lie in i) the selection criteria of member 
galaxies of past studies, which were often not based on exten¬ 
sive spectroscopic information; ii) the dynamical status of the 
analyzed clusters; or iii) an evolutionary trend when considering 
clusters at different redshifts. We defer a more complete analysis 
of this question when our analysis has been extended to other 
CLASH-VLT clusters. The specific case of MACS 1206 is fur¬ 
ther analyzed at the end of this section. 

The spatial separation between PAS and ELG galaxies is also 
clear in Eig. Q (top panel), where we show the distribution of 
galaxies of different spectral classes in the plot of rest-frame 
velocity vs projected clustercentric distance, the so-called pro¬ 
jected phase-space diagram. In the same figure we also plot the 
integral and differential (V) and cry profiles (middle and bottom 
panels). The integral {V) profile shows that the measure of the 
mean velocity is independent of the limiting radius of the sam¬ 
pled region in agreement with a scenario of a unimodal, relaxed 
cluster. As for the differential (V) profiles, the most noticeable 
feature is the small value of (V) of the ELG and TOT populations 
at R'S 2 Mpc. This feature is likely related to the presence of 
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a low-velocity substructure detected at $ R 200 in the NE quadrant 
(see Sect. lO) . 

The integral and differential cry profiles show that the ve¬ 
locity dispersion decreases with R. This trend is the result of the 
cluster mass density profile and the anisotropy velocity profile as 
discussed in B13. The interesting feature is that the face values 
of (Tv as computed for the ELG population are higher than those 
of the PAS population at comparable radii. This is also shown 
in Eig. 3 of B13, where the two populations are color defined. 
However, the errors in individual bins are large and the differ¬ 
ence is scarcely significant. To further investigate this point we 
used two approaches. We fitted the velocity dispersion profiles 
of ELG and PAS galaxies as given in Pig. [T] (bottom panel) ob¬ 
taining lg(cry_ELG) = 3.15(+0.124) - 0.399(+0.135) x lg{R) and 
lg(o‘v,PAs) = 2.98(+0.121) - 0.161(+0.047) x lg(R), where cry 
and R are in units of km s ' and Ii^^q Mpc. The slopes are differ¬ 
ent at the 92% c.l. according to the Welch test (e.g., Guest [T9MT) . 
In the second approach we considered the region between 1 and 
3 Mpc, where both populations are well represented, obtain¬ 
ing cryELG,i -3 = for 99 ELG galaxies and crypAs,i -3 = 

SST+gg for 82 PAS galaxies, i.e., cryELG,i -3 > crypAs,i -3 at the 
97% c.l. according to the P-test (see, e.g.. Press et al. 119921) . We 



0 2 4 6 

R [Mpc] 


Fig. 7. Top panel: rest-frame LOS velocity vs projected cluster¬ 
centric distance for galaxies of different spectral classes. The 
cluster center coincides with the position of the BCG and the ver¬ 
tical dashed line indicates the value of f? 200 - The horizontal green 
line indicates the BCG velocity. Middle panel. Integral profile 
of the mean LOS velocity shown through its error bands at the 
68% C.L. The mean velocity at a given projected radius from 
the cluster center is estimated by considering all galaxies within 
that radius. Small black dots, small red circles, and blue circles 
show the differential profiles for TOT, PAS, and ELG galaxies 
(each point is the value computed using 20 galaxies). Bottom 
panel. As in the middle panel, but for the LOS velocity disper¬ 
sion (here each point is the value computed using 30 galaxies). 
The green square indicates the estimate of the X-ray temperature 
Tx reported by Postman et al. (120121 ) using Chandra data within 
the R - 0.714 Mpc radius, indicated by the position of the 
green arrow, and here converted to cry assuming the density- 
energy equipartition between ICM and galaxies, i.e., jSspec = 1 
(see Sect. EB. 


conclude that, in the case of MACS 1206, the PAS and ELG pop¬ 
ulations do differ in their velocity dispersion but this difference 
is hidden by the combination of the trend of the cry profile with 
the strong PAS vs ELG spatial segregation when global values 
are analyzed (see the first part of this section). 


4. Cluster substructure 

We applied a set of tests in the velocity space (ID-tests), in the 
2D space of positions projected on the sky, and in the combined 
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2Dh- 1D space (3D-tests). Since the sensitivity of individual di¬ 
agnostics depends on the relative position and velocity of sub¬ 
structure, no single substructure test is the most sensitive in all 
situations (Pinkney et al. ll996t . 

4.1. Analysis of the velocity distribution 

The velocity distribution was analyzed to search for possible 
deviations from Gaussianity that might provide important sig¬ 
natures of complex internal dynamics. We used the two robust 
shape estimators by Bird & Beers (I1993I) . the asymmetry in¬ 
dex (AI) and the tail index (TI). Following the Indicator test by 
Gebhardt & Beers (I1991I) . we also checked the peculiarity of the 
BCG velocity (Vbcg = 131 849 km s ') in all the samples con¬ 
taining the BCG (namely the TOT, R200, and PAS samples). 

To detect and analyze possible departures from a single-peak 
distribution, we used the ID-DEDICA method, already applied 
in Sect. 12 to determine MACS 1206 membership. DEDICA was 
introduced by Pisani (119931 1. and extended to the multivariate 
case in Pisani (119961 1. The algorithm is described in the orig¬ 
inal studies and references therein. Here we summarize useful 
points (see also Appendix A in Girardi et al. 119961) . DEDICA is 
a method of cluster analysis based on the estimate of the prob¬ 
ability density of a data sample, where the density is estimated 
by using an iterative and adaptive kernel estimator, in this case a 
Gaussian kernel with an optimal choice of the width. It is based 
on the assumption that a cluster corresponds to a local maximum 
in the density of points. This method allows us to estimate i) the 
statistical significance of each subcluster it detects and ii) the 
probability that each galaxy in the data sample is a member of 
each detected subcluster. In particular, it gives a measure of the 
overlapping between two systems (see Girardi et al. 119961 and 
Pisani et al. 119961 ). We note that DEDICA is a non-parametric 
method in the sense that it does not require any assumption about 
the number of clusters or any other of their features. Table 0] 
summarizes the results of the ID-DEDICA method applied to 
the TOT and other samples of MACS 1206 and Eig.[8]shows the 
relevant plots. 

The velocity distributions of the TOT and R200 samples 
show no departures from a Gaussian, nor signs of a peculiar 
BCG velocity, and only one significant peak is detected by the 
ID-DEDICA reconstruction. Within 0.5R2oo^ two peaks, formed 
of 91 and 37 galaxies, are detected by the ID-DEDICA recon¬ 
struction (see Eig.[8l). However, the two groups are strongly over¬ 
lapped with 26-1-26 galaxies having a non-negligible probability 
to belong to both groups and not different in their 2D distribu¬ 
tions. 

The velocity distributions of galaxy populations per spectral 
class show marginal or no departures from a Gaussian with the 
exception of the HDSr population. In the case of the PAS class 
no sign of a peculiar BCG velocity is found. Moreover, only one 
significant peak is detected by the ID-DEDICA reconstruction 
with the exception of wELG, a very poor class where a non¬ 
significant peak can be detected. The velocity distribution of the 
PAS galaxies outside R 200 and that of the PAS galaxies within 
O. 5 R 200 show the presence of two peaks, but, as in the case dis¬ 
cussed above, the two groups are strongly overlapped and show 
no difference in their spatial distribution. 

We note the peculiar velocity distribution of ELGs within 
R 200 , which is significantly peaked at lower values with re¬ 
spect to the mean velocity of the whole cluster, (Velg,<r20o) = 
(130709+171) km s“'. The velocity distribution of vsELGs is 
peaked around the same value and, in addition, the visual in¬ 
spection of the ID-DEDICA reconstruction (Eig.[8] bottom-right 
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Fig. 8. Galaxy density reconstruction in the LOS velocity space 
through the ID-DEDICA method. Results for all galaxies in the 
TOT, PAS, and ELG samples are shown from left to right, while 
results for the whole and partial cluster regions are shown from 
top to bottom. The result in the vsELG sample is shown in the 
bottom-right panel. In each panel, the dashed vertical line indi¬ 
cates the BCG velocity. Units on the y-axis are rescaled to the 
maximum value of each density distribution. 


panel) suggests a secondary bump at high velocities. These fea¬ 
tures are related to the presence of a low-velocity NE galaxy 
structure and of a high-velocity SEext galaxy structure as we 
discuss in Sect. 14.31 We do not plot the reconstructions for HDSr 
and wELG, whose distributions are poorly sampled. 

4.2. Analysis of the 2D galaxy distribution 

The main features in Eig. 0] are the WNW-ESE elongation and 
the lack of ELGs in the SW quadrant with respect to the NE 
quadrant (14 vs 54 galaxies). In the computation of ellipticity 
(e) and position angle of the major axis (PA), we followed the 
moments of inertia method (Carter & Metcalfe 119801 see also 
Plionis & Basilakos 120021 with weight w - 1). Table |5] lists 
the values of e and PA (measured counterclockwise from north) 
for the TOT and R200 samples, and the main spectral classes. 
Pigure|2 shows the integral estimates of e and PA at increasing 
radii. As far as concerning PAS galaxies, the value of e increases 
out to R ~ R 200 , which is the radius containing most PASs, and 
then is stable. The value of PA is already stable at R^ ^ 200 - As 
for ELGs, Eigure fallows us to appreciate a strong variation of 
e and PA at R ~ 2-3R2oo, while values are roughly stable at R^ 3 
Mpc. This variation is likely due to the NE vs SW quadrant 
asymmetry and, specifically, to the presence of an important NE 
structure as discussed at the end of the section. When comparing 
values at large radii, R^ 3 Mpc, the PAS and ELG spatial 
distributions have similar position angles, PA= 110-120°, while 
the ELG distribution is somewhat rounder than the PAS distribu¬ 
tion, epAs ~ l -5-2eELG at the 2-3cr c.l. depending on the precise 
radius. Table|5]also lists the values of e and PA for the HDSr pop- 
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Table 4. Results of the substructure analysis. 


Sample 


TI 

% 

AI 

% 

ID 

VBcopec. 

% 

DED. 

Np 

2D 

DED. 

Np 

Vgrad. 

% 

2D+ ID 

DS(y) DSo-v,corr (DSo-v) 
% % 

WHOLE 

445 

ns 

ns 

ns 

1 

1 + 3 

ns 

ns 

ns (> 99.9) 

R200 

249 

ns 

ns 

ns 

1 

1 + 1 

ns 

ns 

ns (ns) 

PAS 

183 

90-95 

ns 

ns 

1 

1 + 1 

ns 

ns 

94 (92) 

HDSr 

38 

90-95 

95 - 

99 

1 

2 

ns 

ns 

ns (ns) 

wELG 

17 

ns 

ns 

- 

- 

1 

ns 

- 

- 

ELG 

168 

90-95 

ns 

- 

1 

7 

ns 

> 99.9 

ns (99.7) 

vsELG 

43 

ns 

ns 

- 

1 

2 

ns 

99.1 

ns (94) 


Notes. Column 1: sample ID; Col. 2: number of galaxies, N^; Cols. 3 and 4: significance of the deviations from Gaussian according to the tail and 
asymmetry indices; Col. 5: significance of the peculiarity of the BCG velocity; Col. 6: number of peaks detected through the ID-DEDICA method, 
only peaks detected with a c.l. larger than 99% are considered; Col. 7: number of peaks detected through the 2D-DEDICA method, where nl-tn2 
indicates the presence of nl major peaks and n2 minor with very low-density peaks; Col. 8: significance of the existence of a velocity gradient; 
Cols. 9 and 10; significance of the existence of substructure according to the DS<L>-test and the DScry-test, the latter for both the profile-corrected 
and standard versions. Only significance values larger than 90% are reported, while non-significant values are indicated with ns. 



Fig. 9. Integral profiles of ellipticity {upper panel) and position 
angle {lower panel) for the whole galaxy population (solid black 
line) and per spectral class according to the labels in the upper 
panel. The values of e and PA at a given radius R are estimated by 
considering all galaxies within R. Results within 1 Mpc are 
too noisy to be meaningful. The vertical dashed line indicates 
the value of R 200 - 


ulation, but the respective profiles are not shown in Fig.j^for the 
sake of clarity. The value of PAncsr agrees with that of PApas 
within Icr c.l.. The value of enosr is larger than epAs only at a 
$ 2cr c.l., the large value of ennsr being due to the few galaxies 
of the dense ESE peak (see in the following). 

We also analyzed the galaxy spatial distribution through the 
2D-DEDICA analysis (Pisani I1996t . In Table |6l we present 
the full information for the relevant peaks, i.e., those with a 
c.l.> 99%, with a relative density with respect to the main 
peak ps~ 0.20, and with at least 10 assigned galaxies. The rel- 


Table 5. Ellipticity and position angle of the galaxy distribution. 


Sample e PA 
_ 


TOT 

0 31 + 0-^4 

117 

R200 

- 0.06 

102 

PAS 

0.41+0“ 

111 

HDSr 

0.55+fll 

105 

ELG 

0 26+^:^*^ 
^•^^- 0.05 

121 


evant maps are shown in Figs. [TO] and [TT] When analyzing the 
whole cluster, the 2D-DEDICA density reconstruction confirms 
the cluster elongation and shows a main structure peaked on the 
BCG with a few secondary peaks (at ESE, SE, and WNW). The 
SE peak in the internal cluster region (hereafter SEint), with 38 
assigned galaxies, is the only secondary peak detected in the 
R200 sample. 

When considering PAS galaxies, the main feature is the clear 
elongation along the WNW-ESE direction in the plane of the 
sky, with the presence of an external secondary peak at WNW. 
Noticeably, HDSr galaxies trace a very dense peak ~ 2 h:jQ 
Mpc at ESE, denser than that around the cluster center. This 
ESE peak is formed of only six galaxies, but its existence is sup¬ 
ported by our independent analysis in Sect. l4.4l based on galax¬ 
ies of different colors. We detected no peculiarity in the velocity 
field in the ESE region. The 2D-DEDICA analysis of the wELG 
population shows only one peak roughly centered on the cluster 
center. To the contrary, the 2D distribution of ELGs is clearly 
not radially symmetric, with a lack of galaxies in the SW region 
and shows several clumps (see Table 1^. The ELG distribution 
shown in Eig. [TT] allows us to better interpret the variation of e 
and PA in Pig. 0 The NE structure centered at R'S R 200 is the 
cause of the value of PA~ 50° atR~ Raoo- When more external 
ELG structures, somewhat related to the NE-SW direction, are 
included the PA value shows a fast, strong change and the e value 
also increases. According to our analysis in Sect. 13.21 the spa¬ 
tial distribution of vsELGs is different from that of mELGs and 
sELGs. The reason is that vsELGs are concentrated in the NE 
and SEext structures, as shown by the contour maps for vsELGs 
and mELGs-i-sELG galaxies (see Pig.fTSli. 
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Fig. 11. Spatial distributions and isodensity contour maps for the PAS {upper left panel), HDSr (upper right panel), wELG (lower 
left panel), and ELG (lower right panel) populations. Labels refer to the peaks listed in Table |6] The /?200 region is highlighted by 
the circle centered on the BCG. 


The study of Pisani (119961 1 shows that the 2D-DEDICA 
method, based on the centrally symmetric two dimensional 
Gaussian kernel, can reconstruct the density of an elongated 2D 
system even in the case with an axial ratio alb - 0.1, although 
slightly rounder isoplethes are obtained (see his Pig. 2). Pisani 
( 119961 1 also shows that a simulated cluster with ajb - 0.1 and 
strongly contaminated by a population of randomly distributed 
points (100 random points and 100 cluster points) is split into a 
series of clusters strongly overlapping each other along the ma¬ 
jor axis (see his Pig.4 and Table 1). This second case is very 
extreme with respect to MACS 1206 both for the presence of an 
important background and the high ellipticity (cf. ajb = 0.1 and 


ajb - 0.7). However, to leave no room for doubt in the case 
of MACS 1206, we run an additional test. We performed the cir¬ 
cularization of the coordinates in the TOT sample, leaving the 
cluster center fixed and using the ellipticity listed in Table|5] and 
repeated the full 2D-DED1CA analysis. We confirmed the pres¬ 
ence and the relative density of the main WNW and ESE peaks, 
but not the SEint peak. However, the SEint peak and its relative 
density are confirmed in the re-analysis of the R200 sample af¬ 
ter the circularization procedure. The WNW and ESE peaks are 
also independently confirmed in the re-analysis of the PAS and 
HDSr samples. 
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Table 6 . Results of the 2D-DEDICA analysis. 



Subclump 

As 

a(12000), 5(12000) 
h : m : s, ° : ' : " 

Ps 

As 

TOT - main 

219 

12 06 11.5 -08 47 57 

1.00 

103 

TOT - WNW 

102 

12 05 49.4 - 08 46 44 

0.20 

29 

TOT - SEint 

38 

12 06 24.9 - 08 49 37 

0.20 

16 

(TOT - ESE) 

41 

12 06 35.0 - 08 49 26 

0.15 

15 

R200 - main 

203 

12 06 11.5 -08 47 57 

1.00 

80 

R200 - SEint 

30 

12 06 24.8 - 08 49 40 

0.20 

17 

PAS - main 

126 

12 06 12.2-08 48 05 

1.00 

41 

PAS - WNW 

44 

12 05 51.0-08 46 30 

0.25 

12 

(HDSr - ESE) 

6 

12 06 37.0 - 08 49 37 

1.00 

6 

HDSr - main 

15 

12 06 11.0-08 47 57 

0.67 

4 

wELG 

15 

12 06 12.0 - 08 47 20 

1.00 

7 


ELG - NE 

39 

12 06 21.8 -08 44 00 

1.00 

14 

ELG - NW 

30 

12 05 57.2 - 08 40 38 

0.85 

11 

ELG - WNW 

24 

12 05 42.0 - 08 46 34 

0.68 

9 

ELG - SE 

13 

12 06 36.0 - 08 52 34 

0.58 

5 

ELG - SEext 

20 

12 06 49.7 - 08 54 19 

0.43 

10 

ELG - NEext 

17 

12 06 31.0-08 38 10 

0.31 

8 

ELG - ESEext 

10 

12 06 54.2 - 08 48 28 

0.31 

5 

vsELG - NE 

24 

12 06 20.5 - 08 44 20 

1.00 

9 

vsELG - SEext 

15 

12 06 49.3 -08 53 13 

0.55 

4 


Notes. Column. 1: subsample/peak ID, labels in parentheses corre¬ 
spond to galaxy peaks that do not completely fulfil the criteria we fixed 
for galaxy number or density, but are listed here becasue they are dis¬ 
cussed in the text; Col. 2: number of assigned member galaxies, As; 
Col. 3: R.A. and Dec of the density peak; Col. 4: relative density with 
respect to the highest peak, psl Col. 5: value of the peak. For each 

population, the number of peaks is summarized in Table|4](Col. 7). 



Fig. 10. Spatial distribution and isodensity contours of the 445 
cluster members (TOT sample). Large cross and circle highlight 
the BCG position and the R 200 region. Small crosses indicate the 
secondary density peaks (see Table|^. 



X [arcmin] 



Fig. 12. Spatial distribution and isodensity contour maps of 
vsELGs {upper panel) and other ELGs (i.e., mELGs-i-sELGs, 
lower panel). The position of labels are those in Eig. fTTI (lower 
right panel) to allow an easy comparison with that figure and 
highlight the concentration of vsELGs in the NE and SEext 
structures. The R 200 region is indicated by the circle centered 
on the BCG. 

4.3. Combining position and veiocity information 

The existence of correlations between positions and velocities 
of cluster galaxies is always a strong footprint of real substruc¬ 
tures. To study the 3D cluster structure we applied two tests. The 
presence of a velocity gradient was searched for by performing 
a multiple linear regression fit to the observed velocities with re¬ 
spect to the galaxy positions in the plane of the sky (e.g., den 
Hartog & Katgert [19961 Girardi et al. 119961 1. The significance is 
based on 1000 Monte Carlo simulated clusters obtained by shuf¬ 
fling galaxy velocities with respect to their positions. We found 
no significant evidence for a velocity gradient. 

The A-statistics devised by Dressier & Schectman (119881 
hereafter DS-test) is a powerful test for 3D substructure, which 
is valid in samples down to 30 member galaxies, at least in the 
case of major mergers (Pinkney et al. 119961) . Eor each galaxy, the 
deviation 5,- is defined as - [(Ann + l)/o'y][((y)ioc - (V))^ + 
(cTyioc - CTy)^], where the subscript loc denotes the local quan¬ 
tities computed over the Ann =10 neighbors of the galaxy, (V) 
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and <Tv are the global quantities, and A is the sum of the |5,| of 
the individual N galaxies of the sample. The significance of A, 
i.e., of the existence of substructure, is checked by running 1000 
Monte Carlo simulations, randomly shuffling the galaxy veloci¬ 
ties. Here, we used two kinematical estimators alternative to the 
6i parameter. We considered separately the contribution of the 
deviation of the local mean velocity from global mean velocity 
Si,v - [(A^nn + X ((^)ioc “ (^)) ^^d the contribution of 

the deviation of the local velocity dispersion from the global ve¬ 
locity dispersion (5/_s = [(A^nn + x (cryioc - cry) (Barrena 

al. 1201 11 see also Girardi et al. 119971 ). The results are listed in 
Table El 

The DS{y)-test returns a strong positive detection of sub¬ 
structure in the case of the ELG population (> 99.9% c.l.), with a 
low-velocity group within /? 2 oo at NE and a high-velocity group 
in the external SE region (see Eig. [T3] for the so-called bubble 
plot). To be more quantitative, we resorted to the technique de¬ 
veloped by Biviano et al. (120021 see also Girardi et al. 120061) . 
We compared the distribution of the 5,_y values of the real ELG 
galaxies to the distribution of the 5,;y values of the galaxies of 
all the 1000 Monte Carlo simulated ELG samples (Eig. [14]). In 
agreement with the result of the DS(y)-test, the two distribu¬ 
tions differ at the > 99.99% c.l. according to the IDKS-test. The 
distribution of the values of real galaxies shows a tail at large 
positive and negative 5,_y values, populated by the galaxies be¬ 
longing to the 3D substructure. In particular, the low-velocity tail 
is more populated than the high-velocity tail in agreement with 
the relative importance of the low-velocity NE substructure in 
Eig. [T3] To quantify the number of galaxies involved in the 3D 
substructure, we assumed a threshold of |d,;y| = 3 and counted 46 
real galaxies with |5,_y| > 3, while five galaxies are predicted by 
simulations; in other words, ~ 41 ELGs are expected to belong 
to the 3D substructure. The galaxies with |(5,_y| > 3 are spatially 
separated in the NE and SE clumps allowing us to estimate the 

rest-frame velocity of these clumps, AErf ^E-550 km s^'and 

Ayrf.sEext ~ +550 km s^'. 

We note that the tails of the ELG 5,,y distribution are pop¬ 
ulated by several vsELGs; 20 vsELGs, i.e., about half of the 
whole vsELG population, have |5,;y| > 2.5. This prompted us to 
apply the DS(y)-test to the ELG subclasses. The test applied to 
the vsELG class returns, in spite of the small number of vsELG 
galaxies, a 99.1% c.l. signal of substructure (see Eig. [B]), while 
no signal is found for mELG and sELG. Erom these results and 
those of the previous section we can conclude that vsELGs are 
strongly concentrated in the NE and SEext substructures and are 
related to the 3D substructure. 

In the context of the above results, the features in the ve¬ 
locity distributions of ELGs and vsELGs plotted in Eig. [8] can 
be explained. The velocity distribution of ELGs within /?200 is 
peaked at low velocity owing to the low-velocity NE substruc¬ 
ture, very important within /?200 (see Fig. [T3I) . In the case of 
vsELGs (whole sample), the high-velocity SEext substructure 
is also important resulting in a velocity distribution suggestive 
of two peaks. 

As for the DScry-test, the MACS 1206 dataset is so large 
in galaxy content and in sampled area that the radial cry pro¬ 
file (Fig. |2l bottom panel; see also Fig. 3 in B13) should be 
taken into account. The application of the standard DScry-test in 
the TOT sample detects a spurious strong signal of substructure 
(> 99.9% in the TOT sample) with a peak of high local veloc¬ 
ity dispersion in the central cluster region and low values of lo¬ 
cal velocity dispersion in external cluster regions (see Fig. [Thl) . 
We overcame this difficulty by correcting the kinematical esti¬ 
mator. Specifically, we consider the deviation of cryjoc from the 



Fig. 13. DS(y) bubble-plot for the ELG population. The larger 
the circle, the larger is the deviation of the local mean velocity 
from the global mean velocity. Blue and heavy red circles show 
where the local value of mean velocity is smaller or larger than 
the global value. Labels indicate the two peaks detected in the 
2D analysis (see Table|6]and Fig.[TT] lower right panel). The /? 2 oo 
region is highlighted by the large black dashed circle centered on 
the BCG. 



Fig. 14. The distribution of the ^,_y values for the ELG popula¬ 
tion. The solid line histogram indicates observed galaxies. The 
dashed histogram indicates the galaxies of simulated samples, 
normalized to the observed galaxy number. The blue vertical 
dashed lines highlight the \6iy\ > 3 regions where most galaxies 
are expected to belong to the 3D substructure. 
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velocity dispersion expected at the corresponding radius R, i.e., 
^i,s,corr “ [(Ainn + ly^^/a-y] X (cTyjoe “ ctv^r). The value of cry,R 
was obtained by fitting the values of the cry-profile in Fig. [T] 
[lg(cryR/km s“')= 3.027 - 0.224xlg(f?//!“Q Mpc)]. This correc¬ 
tion cancels the spurious signal of substructure in the TOT and 
other samples (cf. DScrycon and DScry values in Tableland cf. 
upper and lower panels in Fig. [Tbl l. We note that no significant 
signal is obtained in the R200 sample, which is by definition less 
spatially extended, even before the correction. 

As for the PAS sample, the DScry-test substructure signal is 
not radially symmetric (see Fig. [Ill) and, as a consequence, does 
not vanish when applying the correction based on the cry profile. 
However, the evidence of substructure is always marginal (at the 
92-94% c.L, see Tabled and few galaxies are related to this 3D 
substructure. We do not discuss this small feature further. 

All the above results are obtained using the classical version 
of DS-test local groups with Ann = 10 neighbors. We also ap¬ 
plied the version with Ann = 20 for the samples having many 
members (TOT, R200, PAS, and ELG) with no relevant differ¬ 
ence in the results. 

4.4. Substructure analysis as a function of galaxy colors 

We present here the main results of our substructure analysis of 
MACS 1206 by binning the sample on the basis of galaxy col¬ 
ors. We divided the sample into four classes. Starting from the 
red sequence galaxies, i.e., those with |(Z? - Rc)difrl ^ 0.31 0, 
we defined the red sequence galaxies with positive and negative 
(B - Rc)diff, called the RedU sample (72 galaxies) and the RedD 
sample (133 galaxies), respectively. We then defined a Blue sam¬ 
ple, containing 152 galaxies with (B - Rc)diff < -0.5, and at in¬ 
termediate colors the Green sample, including 74 galaxies with 
-0.5 < (B - Rc)diff < -0.31. By repeating the aforementioned 
analysis per color class, we found no strong evidence of non- 


We adopt a threshold of 0.31 instead of 0.3 because we chose to 
include the BCG in the red sequence. 



10 5 0 -5 -10 

X [arcmin] 



X [arcmin] 

Fig. 16. DScTy bubble-plot in the standard and profile corrected 
cases for the TOT sample (upper and lower panels). The larger 
the circle, the larger is the deviation of the local velocity dis¬ 
persion from the global velocity dispersion. Blue and heavy red 
circles show where the local value is smaller or larger than the 
global value. In the upper panel the signal of high local cry in the 
central region and low local cry in the external region reflect the 
trend of the cry profile (see Fig. [7] bottom panel) and is no longer 
significant when applying the suitable correction (lower panel). 
The R 200 region is indicated by the dashed circle centered on the 
BCG. 


Gaussianity in the velocity distribution (i.e., always < 95% c.L), 
and no evidence of multimodal distributions according to the 
ID-DEDICA method. We show the results of the 2D analysis 
in Eig. [18] The 2D density maps for Green and Blue galaxies 
resemble those of wELGs and ELGs in Eig. M RedU galaxies 
are strongly clustered around the cluster center; the WNW peak 
has a density ps well under the threshold we fixed as being rel¬ 
evant and it is displayed for completeness. RedD galaxies trace 
the elongated structure of the cluster and, interestingly, confirm 
the presence of the ESE peak (with 14 objects). 
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Fig. 17. DScTy bubble-plot as in Fig. [16] (standard case), but for 
the PAS population. 






Fig. 18. Spatial distributions and isodensity contour maps for the 
four color classes; RedU {upper left panel), RedD {upper right 
panel). Green {lower left panel), and Blue {lower right panel) 
populations. Labels are the same as in Fig.fTTI 


5. General results and discussion 

5.1. Cluster structure 

The statistical results of the substructure tests applied to the TOT 
and R200 samples indicate that MACS 1206 does not show the 
presence of large-scale substructure. 

The 2D elongation is the most important large-scale feature 
in the cluster structure. In the plane of the sky the cluster struc¬ 
ture appears clearly elongated along the WNW-ESE direction 
(see Eig. [Toll. We obtained PA = 102);g and ll?);^ degrees for 
galaxies in the R200 and TOT samples, respectively. This ori¬ 
entation matches that of the mass distribution from strong grav¬ 
itational lensing (see Pig. 1 in Zitrin et al. 120121) . BCG shape. 


X-ray isophotes, and the mass distribution from weak gravita¬ 
tional lensing (PA in the range 105-119°, Table 6 in Umetsu et 
al. 120121) . the asymmetric distribution of the intracluster light 
(PA of BCG and ICE in the range 101-109°, Presotto et al. l20T?l l. 
Our novel result is that there is neither evidence for a velocity 
gradient along the WNW-ESE direction nor for a velocity differ¬ 
ence between the main system and most subclumps. This sug¬ 
gests that the main MACS 1206 structure is lying in the plane of 
the sky. 

As for the interpretation, it is well known that cluster’s dy¬ 
namical activity is strongly correlated with the tendency of clus¬ 
ters to be elongated and aligned with the nearby clusters (Plionis 
& Basilakos l2002t because of the anisotropic merging along the 
LSS filaments as also shown by the analysis of simulations (e.g., 
Basilakos et al. 120061 Cohen et al. 120141 ). However, it is also 
known that the final distribution after violent relaxation may 
be aspherical (Aaarseth & Binnev [T978l White [199^ . Thus, in 
principle, the elongation of MACS 1206 is not in contrast with 
the lack of statistical evidence for the presence of important sub¬ 
structure. 

An interesting global parameter is the value of the ratio be¬ 
tween the energy per unit mass of galaxies to that of ICM as 
parametrized with jSspec = cryl{kTxlpmp), where Tx is the X- 
ray temperature, p = 0.58 is the mean molecular weight and 
rup the proton mass. The value ySspec = 1 indicates the density- 
energy equipartition between ICM and galaxies. The mean ySspec 
value observed for massive clusters is consistent with unity both 
in nearby (Girardi et al. 119961119981 ) and distant systems out 
to z ~ 0.4 (Mushotzky & Scharf 120041 ). while large devia¬ 
tions may be related to important cluster merging phenomena 
(e.g., Ishizaka 119971 ). In the case of MACS1206, the value of 
kTx - 10.8+0.6 keV has been obtained by Postman et al. (120121 ) 
using Chandra data within a radius of 0.714 Mpc and ex¬ 
cluding the very central cluster region. Using our estimate of cry 
within the same radius, we obtainedjSspec = 0.96+0.14, in agree¬ 
ment with the scenario of density-energy equipartition (see also 
Fig. [21 bottom panel). 

At small scale, the analysis of the 2D galaxy distribution 
shows three secondary peaks with low values of relative den¬ 
sity ps~ 0.2 with respect to the main cluster peak (WNW, SEint, 
ESE in the TOT sample, see Table[6l). The ESE and WNW peaks 
have already been detected as significant overdensities by Lemze 
et al. (120131 ). who have used a different sample and method. In 
the R200 sample, the SEint density peak is the only one detected, 
leading to a 12% of cluster members (8% in luminosity) related 
to the substructure. Among clusters with no major substructure, 
the typical fraction of galaxies in substructure is ~ 10% (e.g., 
Girardi et al. 119971 Guennou et al. 120141 ). In the TOT sample, 
which extends out to R ~ 2 R 200 we formally assigned 52% of 
the galaxies to substructure. However, we stress that the pre¬ 
cise number of objects assigned to the WNW clump is largely 
uncertain owing to possible border effects. The WNW clump 
contains galaxies out to the western boundary of the sampled 
region. We checked this point extending the size of our sample 
out to ~ 4 R 200 with the addition of a random background field. 
The addition of a field equal to 5% of the TOT sample does not 
change the value of the relative density ps.wNW of the WNW 
peak, but reduces the number of the assigned WNW members. 
The fraction of WNW members in the TOT sample changes from 
WNW/TOT~ 23% to WNW/TOT~ 10% leading to an impor¬ 
tant reduction in the estimated fraction of galaxies in substruc¬ 
tures. In any case, our result agrees with substructure being more 
important in the external cluster regions (e.g.. West & Bothun 

[ 1223 . 
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ESE* 


NW 


Fig. 19. Subaru Rq image of MACS 1206 with superimposed features summarizing the results of the cluster structure analysis. 
Circles, in increasing size order, indicate the density peaks detected in the 2D galaxy distribution of TOT (cyan), PAS (red), HDSr 
(magenta), wELG (green), ELG and vsELG (blue) samples (see Table|^and Eig.fTTTi: the NE and SEext clumps being also detected 
through our 3D analysis. The lines indicate the likely directions of cluster accretion. Blue contours are the smoothed contour levels 
of the X-ray surface brightness derived from Chandra archival data (in the 0.5-2 keV energy range; see also Ebeling et al. 12009b . 
The large black circle encloses the f ?200 = 2 Mpc ~ 5.9' region. 


5.2. Tracing a past phase of accretion 


Both the PAS and HDSr populations trace a cluster structure 
elongated along the WNW-ESE direction suggesting that this is 
the direction of the main accretion filament during the cluster 
formation. Our interpretation is that HDSr galaxies represent the 
galaxies more affected by the last important phase of accretion, 
likely a past merger or mergers along the WNW-ESE direction. 
In fact, the HDSr spectral features reflect the evolution of galax¬ 
ies having undergone starburst or star-forming activity, which 
has been suppressed by some physical process, possibly related 
to the high-density environment, and observed 1-2 Gyr later (see, 
e.g., M04). This time is comparable with the crossing time of the 
present cluster (R 2 oo/rrv ~ 1.8 Gyr) indicating that the supposed 
merger(s) are in an advanced phase. This also agrees with the 
fact that we do not see relevant features of this past dynami¬ 
cal activity, since the populations of the intervening subclusters 
have had time to mix. In principle, we cannot exclude a phase of 
violent relaxation, since its time scale is comparable to the dy¬ 
namical time. A time estimate of 1-2 Gyr is also consistent with 
the elongation of the X-ray isophotes which is a long-lived phe¬ 
nomenon (e.g., Roettiger et al. 11996b and the asymmetry of the 


ICE, whose surviving time can be comparable with one crossing 
time (Rudick et al. 12009b . 

The observational evidence in favor of a past merger along 
the WNW-ESE direction is the presence of a secondary dense 
peak in the distribution of HDSr galaxies (ESE peak), conhrmed 
in the distribution of RedD galaxies, too (see Eig.[TT]and Eig.fTSl 
upper right panels). The mean values of EW(H(5) of the RedU 
and RedD populations (0.19 + 0.35 A and 1.10 + 0.24 A) are 
only 2cr different, but in the direction of reinforcing our idea that 
the ESE peak is only traced by the younger population (RedD). 
The HDSr ESE peak might be the remnant of a galaxy system, 
now merged to form MACS 1206. However, it is formed only of 
faint galaxies (Rq >21.5 mag), lacking a dominant galaxy which 
would have supported its identification with the remnant core of 
an important subcluster. Another point in favor of a past merger 
along the WNW-ESE direction is the presence of a bright galaxy 
[at R.A.^12'»06‘ "15?66 , Dec.=-08° 48'21.8" (J2000); also G2 
in Grillo et al. I2014bll . which is the second brightest galaxy 
within 0.5 Mpc from the cluster center, and which lies at 
the extremity of the asymmetric ICE distribution, and seems to 
have strongly suffered from the interaction with the BCG and/or 
the central cluster potential (see Eig. 6 of Presotto et al. 12014b . 
This galaxy is also just at the border of our inclusion in the HDSr 
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sample [EW(H5)=3 + 1 A]. A third point is the WNW overden¬ 
sity in the spatial distribution of PASs, although the presence of 
a close overdensity of ELGs makes the interpretation difficult. 

To understand the physical mechanisms generating HDSr 
galaxies, we note that, at large scale, their spatial distribution 
follows that of PAS galaxies. This agrees with the findings of 
Dressier et al. (I2013I I. who also suggested an evolutionary path 
of the type PAS=>PSB=>PAS. However, we cannot completely 
agree with their scenario since our morphological analysis based 
on Sersic indices shows that HDSr galaxies are intermediate be¬ 
tween PAS and ELG galaxies (see Eig.[3]l. The intermediate mor¬ 
phology suggests that HDSr galaxies are the result of two mixed 
populations or that the physical process involved in their trans¬ 
formation works on both their gas content and morphology. The 
second interpretation favors a scenario of merger or tidal inter¬ 
action rather than a gas shocking or stripping. In particular, the 
dense HDSr ESE concentration could be explained by a mech¬ 
anism similar to that suggested by Struck (120061 1. According 
to this mechanism, the gravitational pull of the cluster core on 
a small group falling through the core will cause the group to 
shrink, increasing its density by an order of magnitude, and as a 
result increasing the galaxy merger rate. 

5.3. Recent and ongoing galaxy infall 

Thanks to the large spatial coverage of the spectroscopic sample, 
the ELG population can be studied in detail. We identify several 
clumps (see Eig. (TT] right lower panel). These galaxy clumps 
are very loose and should be interpreted as overdensities in the 
LSS filaments accreted to the cluster rather than gravitationally 
bound, isolated groups. These structures are generally external 
to R 200 with the exception of the NE clump which, being charac¬ 
terized by a AVNE.if-550 km s“', might be seen there because 

of projection effects. The ELG clumps we detect are likely the 
cluster building blocks during the ongoing galaxy infall. 

N-body simulations by Balogh et al. (120001 1 and Moore et al. 
( 120041 1 follow the so-called backsplash or overshoot of galaxies 
that have passed through the cluster and joined the virialized sys¬ 
tems: both these studies identify the zone between 1 and 2 R 200 
as the overlap region where roughly half the galaxies are mem¬ 
bers of the cluster and half are infalling. Since individual ELG 
galaxy clumps are still detectable, we expect that their mem¬ 
bers are likely infalling rather than backsplashing. Considering 
the region between 1 and 2 R 200 , the ratio between galaxies 
that belong to the ELG structures and the whole population is 
93/161 ~ 0.6. Our result suggests that 60% of galaxies have 
never passed through the cluster in rough agreement with the 
50% predicted by simulations by Balogh et al. (120001 1 and Moore 
et al. (|2004] |. 

When analyzing the geometry of the infall, the ELG clumps 
provide a few preferred directions, traced in Eig. [19] roughly 
connecting the clumps and the cluster center. Figure [T^ shows 
how both the WNW-ESE and NW-SE directions are somewhat 
matched in the elongation of X-ray isophotes, which are tilted 
going from the internal to the external regions. The net result of 
the secondary infall is a somewhat rounder assembling than that 
due to the past WNW-ESE accretion as suggested by the differ¬ 
ent ellipticity between PAS and ELG galaxies (see Fig. 0 and 
Tablejg. 

The galaxies of the vsELG class are spatially segregated with 
respect to those of mELG-HsELG classes, clustering in the NE 
and SEext clumps. Figure |20| shows the values of stellar mass 
and SFR that Annunziatella et al. (I2014I I computed for cluster 
member galaxies. Comparing vsELGs and other ELGs in the 



log(M./M.J 

Fig. 20. SFR vs stellar mass of cluster galaxies. Individual val¬ 
ues are those computed by Annunziatella et al. (I2014I I and are 
shown here for different spectral classes. Black solid circles and 
small squares indicate biweight means with bootstrap errors for 
vsELGs and mELGs-i-sELGs considering three bins per order of 
magnitude in stellar mass. 

10 ®< < 10 '° range, we note that the specific SFR of 

vsELGs is higher, particularly at low stellar masses. We inter¬ 
pret the observed spatial segregation as the result of the same 
physical process simultaneously acting in most galaxies within 
the parent clump, which is enhancing SF in the NE and SEext 
clumps, or, alternatively, quenching SF in other clumps. In order 
to understand whether this process is connected to the cluster or 
rather to a pre-infalling environment, we analyzed the relative 
fraction of vsELGs in the cluster and in the field. To this end, we 
selected a representative set of 776 field galaxies as non-member 
galaxies in the 0.2 < z < 0.7 redshift range, with spectral clas¬ 
sification. We compared the fraction of vsELGs with respect to 
all ELGs. We obtained that the vsELG/ELG ratio in the cluster, 
43/168 (26%), is completely consistent with the one in the field, 
160/602 (27%). This does not support a scenario where cluster 
environment is relevant in determining the vsELG/ELG ratio. 

In the context of possible infalling galaxy structures, we 
stress that our study focuses on galaxies belonging to the main 
cluster peak in the redshift distribution, and does not consider 
foreground or background redshifts peaks. In the study of Young 
et al. (I2014I I. a few galaxies in a foreground redshift peak are re¬ 
lated to the presence of a significant interesting SZE feature, and 
interpreted as a foreground or infalling group. 

5.4. Intermediate galaxy populations 

Here we discuss the two intermediate populations, wELG and 
HDSb, which according to our interpretation represent the fate 
of recently in-fallen galaxies. The HDSb galaxies in our study 
are comparable to the PSBs studied by Muzzin et al. (120141 1. 
which have very strong EW(H5) [EW(H5)> 5A for the average 
stacked spectrum] and are defined with D„(4000) < 1.45, i.e.. 
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Fig. 21. Projected phase space as in Fig. [T] top panel, but with 
normalized units in position and velocity axes. Small magenta 
crosses and large blue crosses indicate HDSr and HDSb galax¬ 
ies. Four of the five HDSb galaxies are located within or close 
to the same ring-like region traced by the cluster PSBs analyzed 
by Muzzin et al. (120141 1: see text for details. 

are blue PSBs (see Fig. 2 of M04). Figure |2T] shows that four 
of the five HDSb galaxies in our sample are located within or 
close to the same ring-like region in the projected phase space 
indicated by Muzzin et al. (120141 see their Fig. 1). In that study 
PSBs are explained with a rapid quenching of the order of $ 0.5 
Gyr after making their first passage of R ~ Q.5R2oo- The same 
time scale of ~ 0.5 Gyr is invoked by M04 after the interrup¬ 
tion of the SF in the HDSb galaxies. In MACS 1206, the five 
HDSb galaxies are elongated towards the SE direction, approx¬ 
imately in the NW-SE direction of the putative accretion flow 
(cf. Eig. 0] and Eig. fTOl l. The wELG population appears distinct 
from the ELG population, since wELGs populate inner cluster 
regions, and from the PAS and HDSr populations since wELGs 
do not follow the WNW-ESE elongation. Our interpretation is 
that wELGs are associated with a more advanced phase in the 
ELG infalling process. 

6. Summary, conclusions, and perspectives 

To elucidate the scenario and physical mechanism leading to 
the assembly of galaxy clusters, and specifically the transfor¬ 
mation of star-forming held galaxies to passive cluster galaxies, 
we present a detailed spectrophotometric study of the massive 
cluster MACS 1206 at z = 0.44 as part of the CLASH survey. 
Our analysis is based on an extensive spectroscopic dataset of 
1920 galaxies over a held of 640 arcmin^ around MACS 1206, 
which was the first observing campaign conducted with VIMOS 
at the VLT as part of the “CLASH-VLT” ESO Large Programme 
including 13 CLASH clusters in the southern sky. The dataset 
is supplemented with multiband photometry based on Subaru 
Suprime-Cam high-quality imaging. We combine galaxy veloc¬ 
ities and positions to select 445 cluster galaxies based on kine¬ 


matic information alone. Using spectral absorption and emission 
features, as well as galaxy colors, we classify 412 galaxies in 
seven spectral classes ranging from purely passive galaxies to 
those with strong Hd absorption (red and blue), and emission 
line galaxies (four classes from weak to very strong lines). We 
then analyze the cluster structure and distribution of these differ¬ 
ent galaxy populations in projected phase-space out to a radius 
~ 2/?200 (4 Mpc). Our specific results can be summarized as 
follows: 

- MACS 1206 appears as a large-scale relaxed system, partic¬ 
ularly within R 200 , with a few, low overdensities detected in 
the projected galaxy distribution (relative density $ 20%). 
The whole galaxy population does not show velocity gradi¬ 
ents or evidence of significant deviations in local mean ve¬ 
locities. The main large-scale feature is the elongation along 
the WNW-ESE direction (PA -117°). 

- Passive galaxies and red, strong Hd galaxies trace the clus¬ 
ter center and the WNW-ESE elongated structure. The lat¬ 
ter also cluster around a secondary, very dense peak 2 
Mpc at ESE. 

- The emission line galaxies are mostly located in several 
loose structures, mostly outside /? 200 - Two of these struc¬ 
tures are also detected in the 3D projected phase space hav¬ 
ing LOS rest-frame velocities AUrf^E ~ -550 km s^’and 
AVrf.sEext ~ +550 km s^'. 

- The clumps of emission line galaxies trace several directions, 
with the net result of a more symmetrical distribution than 
that traced by early-type galaxies (the difference in ellipticity 
is a factor of 1.5-2). 

By studying the spatial and velocity distributions of galaxy 
populations of different spectral types, we can piece together a 
likely scenario for the galaxy assembly history of this massive 
system. The oldest and main component of MACS 1206, made of 
early-type galaxies, clearly traces the main phase of cluster for¬ 
mation. However, a more recent anisotropic accretion phase that 
occurred 1-2 Gyr before the epoch of observation is apparent in 
the distribution of red strong Hd galaxies along the prominent 
WNW-ESE filament, including a strong overdensity of these 
evolved post-SF galaxies (ESE peak), possibly the remnant of 
an accreted group. In addition, we identified a secondary promi¬ 
nent young component in the external regions of MACS 1206 
populated by late-type galaxies, which trace the ongoing infall. 

Our study shows that large spectroscopic samples of cluster 
galaxies extending to ^ 2 R 200 allow us to link the cluster sub¬ 
structure with distinct galaxy populations, whose spectral types 
can be used to age-date their underlying stellar populations and 
SE history. This ultimately allows us to reconstruct the evolution 
of galaxies and their parent cluster, thus going beyond the sim¬ 
plistic concept of radial accretion and following cluster forma¬ 
tion through their building blocks. We plan to extend this anal¬ 
ysis to other clusters in the CLASH-VLT survey, which include 
systems of different dynamical states and over an interesting red- 
shift range. A more refined spectral classification will also be 
possible in clusters in this program, which include -7000 spec¬ 
troscopic members. 
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